in 



c^ 



Slow motion of particle systems as a limit 

of a reaction-diffusion equation with 

half-Laplacian in dimension one 



_ M.d.M. Gonzalez * 

(N ■ R. Monneau t 



July 6, 2010 



< 

. . Abstract 

"^ , We consider a reaction-diffusion equation with a half-Laplacian. In the case where the solution 

is independent on time, the model reduces to the Peierls-Nabarro model describing dislocations as 
transition layers in a phase field setting. We introduce a suitable rescaling of the evolution equation, 
using a small parameter e. As e goes to zero, we show that the limit dynamics is characterized by 
^ . a system of ODEs describing the motion of particles with two-body interactions. The interaction 

^D , forces are in 1/x and correspond to the well-known interaction between dislocations. 
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t^^ ■ Keywords: non-local Allen-Cahn equation, reaction-diffusion equation, singular limit, Peierls- 

p2 . Nabarro model, transition layer, particle systems. 



^ : 1 Introduction 



Let us set the one-dimensional half-Laplacian operator L defined on functions w G L°°(R) fl 
C;jf (M) for some < /3 < 1 by the Levy-Khintchine formula 

Lw{x) = - — {w{x + z) - w{x) - Zw'{x)l{\z\<i]) . 



Given a suitable l-periodic potential W which satisfies in particular W'{k) = for /c G Z 
{W will be precised later in assumption (A)), it is known (see Section [3l2l) that there exists 
a unique solution </> of 

Lcj) - W'{(l)) =0 on M, 
' ' ' "^ ^' > and 0(-oo) = 0, 0(0) = -, 0(+cx)) = 1. 
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Model ( II. ip is known as the Peierls-Nabarro model for dislocations in crystals where W is 
called the stacking fault energy. Here 0(a;) is a phase transition representation of a dislocation 
localized at the position x = 0. For a physical introduction to the model, see for instance 
[25]; for a recent reference, see [37]; we also refer the reader to the paper of Nabarro [33] 
which presents an historical tour on the model. 

We now consider a scalar function v solution of the following evolution problem associated 
to the Peierls-Nabarro model: 

(1.2) vt = Lv-W'{v) + a%t,x) on (0, +oo) x M 

where o"^ is the exterior stress acting in the material. This equation has been for instance 

considered in [32] (see also [T^ for a similar model). We assume that this exterior stress is 

given for e > small, by 

f / \ f t x^ 

a {t,x) =ea\—,- 

where u is a given suitable function. Then we can consider the following rescaling 

f/ N f t X 

v\t,x) =v[—,- 

We now look for scalar functions v'^ satisfying 

vl = -iLv' - -W'{v') + a(t, x) \ on (0, +oo) x 
^1.3) ' ^ ' ^ 

w^(0,x) = Vq{x) for X G M. 

When v^ has A^ transition layers as e goes to zero, we will see that we can associate a 
particle to each transition layer such that the dynamics of v'^ is driven by the following ODE 
system for particles (a;j(t))j=i^...^Ar: 



^1.4) { dt ^ 

Xi(0) = x. 

where 



-o-(t,Xi) + -y^ , on (0, +CX)) 

„o 



for i = l,...,N 
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1.1 Statement of the main result 

To state precisely our main result, we make the following assumptions: 

(A) 

i) Assumption on the potential 



' W e C2./3(M) for some < /3 < 1, 
W{v + 1) = W{v) for any w G M, 
i W = on Z, 



VT > on M\Z, 

^ W"{0) > 0. 
ii) Assumption on the exterior stress 

a E BUC ([0, +oo) X R) and for some < /3 < 1 and K > Q, we have 

|o"z|l°°([0,+oo)xM) < -f^, |o"i|L°o([o^+oo)xR) < -f^, 

\ax{t,x + h) - ax{t,x)\< K\h\l^ for aU x, /i G M, tG[0, +oo). 

We recall that BUC is the space of bounded uniformly continuous functions. 
iii) Assumption on the initial condition 



For x\ < x\ < ... < x%, we set 



:i.6) 



vl{x) 



-a{0,x) + Y,<P\^^^j with a = W"{0)>0. 

i=l ^ ^ 



Let us now give our main result (which has been announced in |15]): 

Theorem 1.1. (Convergence to the ODE system) 

Under assumption (A), there exists a unique solution (xj(t))j=i^...^7v of ( [i.^[ j. Let 

N 

v\t,x) = Y,H{x-x,{t)) 

i=l 

where H{x) = l[o,+oo)(3;) is the Heaviside function. Then for any e > 0, there exists a unique 
viscosity solution v'^ of M.^] . Moreover as £ — )■ 0, v^ converges to v^ in the following sense: 



(1.7) 

and 
(1.8) 



limsup v%t',x') <{v°)*{t,x) 



liminf v'(t' ,x') > (v°)Jt,x) 

(t',x')~^(t,x), e^O 



for{t,x) e [0,+oo) X M. 

Remark 1.2. We recall here that the semi- continuous envelopes of a given any function u 
are defined as follows 



u*{t,x) = limsup u{t',x') and u^{t,x) = liminf u{t',x'). 

(f ,x')^(t,x) {f ,x')-^(t,x) 



The proof of Theorem 11.11 is done by constructing suitable sub and supersolutions essen- 
tially based on the following ansatz of the solution: 

(1.9) nt, x) = '-ait, x) + J2[^ I^^Z^^ _ exm ^ {^^—^^ } 

where z/; is a suitable corrector. 

1.2 Brief review of the literature 

As we have mentioned, the diffusion equation (11. 3p is the evolution equation for the Peierls- 
Nabarro model, which describes the dynamics of dislocation defects in crystal. For general 
references, we refer the reader to [3] (see also the paper |36j). The Peierls-Nabarro model 
can also be derived from Frenkel-Kontorova models (see |19]). 

In the time-independent setting, our model is related to a pioneering work of Alberti- 
Boucchitte-Seppecher. They have shown in [T| that if we consider the energy functional for 
VL C M" and a double well potential W given by 



Ee[u] = ]: f \Vuf + - [ W{U) 

"^ Jn ^ Jan 



then, in particular, the rescalings E^[u\/ |log£:| Gamma-converge when e — )> to -7/''(S'„q) if 
n = 2 (thus dQ is one dimensional) and the minimizers u^ converge to Uq, where Mq = 0, 1 
on dQ, Auq = in fi, and the singular set is 

Sua = {z & dQ : Uq jumps at z} . 

A generalization to the case of elasticity equations with application to dislocations has been 
done by Garroni and Miiller in 



For the time- dependent case, the (anisotropic) mean curvature motion has been obtained 
at the limit (in the framework of viscosity solutions) by Imbert and Souganidis |26) . Related 
to this result, let us also mention the work |12| . 

It is important to mention that a result similar to Theorem 11.11 has been obtained for 
the non linear heat equation (after a suitable rescaling) 

vt = e^v^^ - W'{v) 

Here the interaction force between particles is also related to the behavior of the layer solution 
associated to the non linear heat equation, and hence is of exponential type. For such results 
using an invariant manifold approach, we refer the reader to Carr-Pego |10| . Fusco-Hale |21] . 
and Ei [16] for generalization to systems of PDEs. For results using the energy approach, see 
Bronsard-Kohn |6], Kalies, Van der Vorst, Wanner [27]. Remark that our approach by sub 
and super solutions seems new (even in this context). We also refer to the paper of Chen 
|llj . for an interesting tour of the different regimes of the solution to the non linear heat 
equation that appear for general initial data. 

Let us mention that similar results have also been obtained for Cahn-Hilliard equations 
or their generalization to systems, called Cahn-Morral systems (see Grant |24]). 

4 



We expect (even if we have no proof) Theorem 11.11 to be true for any power of the 
Laplacian (—A)", a G (0,1), once the layer solutions of (II. ip are found. The existence of 
such layer solutions is the content of the work in preparation [7], that tries to generalize the 
work of Cabre and Sola- Morales J3 for layer solutions of the half Laplacian. The extension 
problem (see |9]) that would be needed in this case is a degenerate elliptic equation with A^ 
weight, that has been well understood in the classical reference |18] . 

We should mention the result by Forcadel-Imbert-Monneau [20], where they homogenize 
in particular system (11.41) and more generally non-local first order Hamilton- Jacobi equations 
describing the dynamics of dislocation lines. Similarly, homogenization results have been 
obtained in [31] for equation (II. 2p for periodic stress a. 

1.3 Organization of the paper 

In Section 2 we review of the notion and properties of viscosity solutions for this type of 
non-local equations. Section 3 contains all the preliminary results that will be required later 
(layer solutions, corrector) together with some motivation of our result and a study of the 
ODE system of particles (II. 4p . Section 4 is the proof of the main result (convergence is 
proven by finding suitable sub and supersolutions of the PDE). Section 5 and Section 6 are 
almost independent of the rest of the paper: Section 5 is a necessary review of known and 
further results (used in Section 6) for the half-Laplacian operator, while Section 6 contains 
the proof of the properties of the layer solution (Theorem 13. ip and of the existence of the 
corrector (Theorem 13. 2p . 

2 Viscosity solutions 

We recall the definition of a viscosity solution v for equation (II. 3p with e = 1 to simplify, 
namely for T e (0, -|-oo): 

{vt = Lv-W\v)^cy{t,x) on (0,T)xM 
[^ 1^(0, x) = Vi^{x) for X G M. 

Definition 2.1. (Viscosity sub and supersolutions) 

An upper (respectively lower) semicontinuous function v G L°^([0,T) x M) is said to be a 
subsolution (resp. a supersolution) of Ii2.10\) if and only if 

v{0, ■) < V* on R 

{resp. v{0, ■) > Vq on M) 
and for any C^ test function ip satisfying 

V = ip at the point (to,a^o) ^ (0,^) ^ ^^ 

and 

V < p on a neighborhood of {to} x R 

{resp. V > p on a neighborhood of {to} x ^) 

then we have 

Pt{to,xo) < L^[p{to,-),v{to,-)]{xo) -W'{p{to,xo)) + a{to,xo) 



where for any functions $(a;), w{x) we set 

L^[^,w]{xo) = - -^mxo + z)-^{xo)-z^'{xo)) + - -^ {w{xo + z) - w{xo)) 

TT J[-l,l] Z-" 71 Jk\[-1,1] ^ 

A function v G L°°([0,T) x R) is said to he a viscosity solution of Ii2.10\) . if and only if v* is 
a subsolution and v^ is a supers olution. 

Finally, a function v G L^^([0, +oo) x M) is said to he a viscosity solution on [0, +oo) x M, 
if and only if its is a viscosity solution on [0, T) xM. for any T > 0. 

Remark 2.2. (Classical solutions are viscosity solutions) 

From this definition, we can easily check that every function v G L°^([0, T) xR)nC;Q'f ([0, T) x 
M) for some < /3 < 1, is indeed a viscosity solution if and only if it is a classical solution 
of li2.1U\) . This is mainly due to the fact that Vt and Lv{t, ■) are then defined pointwise. 

Without entering in the details of the viscosity solutions (for which we refer the interested 
reader the paper by Barles-Imbert |1]), let us recall the known results that we will use. 

Proposition 2.3. (Perron's method) 

Assume that there exists a supersolution v and a subsolution v of equation 1^2. 10\) satisfying 

v<v on [0, T) X M. 
Then there exists a viscosity solution v of Ii2.10\) satisfying 

v<v<v on [0, T) X R. 
Let us also recall the comparison principle: 

Proposition 2.4. (Comparison principle) 

Assume that v is a subsolution (resp. v is a supersolution) of equation ^2. 1 U) on [0, T) such 
that 

t;(0,x) < tJ(0,x) for all x G R. 

Then 

v<v on [0,T)xR. 

In particular, when the initial data Vq is continuous, the comparison principle implies the 
uniqueness and the continuity of the solution. 

3 Preliminary results 

In this section, we collect several results that will be used in the section H] for the proof of 
our main result. 



3.1 Formal ansatz 

In the following, we try to explain formally why we have to expect an ansatz as in (11 .9^ for 
the solutions of fll.3p . 

Namely, we are looking for an ansatz for the solutions of 
(3.11) -evi = -{W'{v')-eLv' -ea} on (0, +cx)) x R. 

Assume to simplify that a is a constant here, and that there is only one transition layer. 

First try for the ansatz 

As a first try, in the case where a is constant, we could consider the Ansatz 

v'{t,x) =ea+-'-' ^^ 



e 
where is the transition layer solution of 

(3.12) L(l)-W'{(j)) = 0. 
When we plug this expression in (13.111) . we get 

XKp' ~ -{W'{ea + (t))-L(j)-ea} 

~ -{W'{ea + (f))-W'{(f))-6a} 

~ W"{(j))a-a + 0{e) 

where we have used (I3.12p for the second line, and a Taylor expansion to get the third line. 
From this computation, we learn two things. First for ^ we have to choose 

(3.13) W"{0)& = a. 

The second information, is that it is impossible to satisfy the equation with this ansatz, and 
we need to introduce a corrector. 

Second try for the ansatz 

Let us take 

-«■/ X ~ , /'x — xi(t)\ , f X — xAt) 

for 5" satisfying (I3.13P and for some constant c to fix later. When we plug this expression in 
(EUD, we get 

Xicj)' + h.o.t. ~ - {W\sa + — ecip) — L(j) + ceLip — ea} 

- - {W'{ea + - ecip) - VT' (0) + ceLi^ - ea} 

~ W"{(j)){a -cip) -a + cLi/j + h.o.t. 
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where we have used a Taylor expansion to get the last line. Finally using (I3.13p . we see that 
ip has to solve 

(3.14) Lij - w"((j)U = —(j)' - -(W'u) - wm. 

c c 

We see that it is convenient to choose (in order to get later ip independent on xi) 

Xi = c. 

Moreover, multiplying 03.14p by 0' we get formally by integration: 

0' (L^ - W"{<P)^) = hct^'f --( [ ^{W'{<P{x)) - W"{0) [ <P' 

./ro c \j]g ax J^ 

i)'f + -W"{0). 
On the other hand using the fact that L is self-adjoint, we get formally 

^'{Lij-W"{^)ij) = f ij{L4>' -W"{^)^') 

Jr 

= 0. 

where the last equality comes from the fact that solves (I3.12p so then 0' solves 

L(j)' - W"{(f))(f)' = 0. 

Therefore, we get 

c = —70" with 7 = 

and tp solves 

L^/j-W"{(f))^lj = (j)' + 7]{W"{(j))-W"{0)) with r] 



'1 



-fW"{0) 



3.2 The layer solution and the corrector 



In this subsection, we state without proofs some results on the layer solution and the cor- 
rector. These results will be proven in Section [6l 

We recall that the layer solution satisfies (II. ip . namely 

r L(j) - W'{(t)) = on M, 
*- ■ '' j 0' > and 0(-cx)) = 0, 0(0) = -, 0(+oo) = 1 

Then we have 



Theorem 3.1. (Existence and properties of the layer solution) 

Under assumption Ai), there exists a unique solution of \3.15\) . Moreover cf) ^ C^ 
some < /3 < 1, and there exists a constant C > such that 

C 
(3.16) < (f)'{x) < —^ — - for all x G M 



l,/3 
loc 



1 + X^ 



and 
(3.17) 



0(x) - H{x) + 



aiTx 



C 

< —- for all \x\ > 1 
x^ 



where H is the Heaviside function and a = W"{0). 

We are also looking for a corrector ip, solution of the following equation 

L'ip-W"{<p)'ip = (l)' + r]{W"{(l))-W"{0)) on R, 



(3.18) 

with 

(3.19) 

Then we have 



i/j{—oo) = = ip{+oo) 



Jm. 



k'\2 



W"{0) ' 



Theorem 3.2. (Existence and properties of the corrector) 

Under assumption Ai), there exists a unique solution tp G H2( 
tp G C^'f (M) n L°°(R) for some < /S < I and 

|V^'|l°°(R) < +00. 



for 



of 1233) ■ Moreover 



(3.20) 



3.3 The ODE system 

Recall that we consider a solution (xj(t))j=i^,,,^Ar of (II ■4p . namely 

^ = ^(-"^"^''^ + iS^)' °" ^°'+^)' for ^ = 1,...,N 
x,(0)=xO, 

where 7 > is given in (II .Sp . 
We recall the following result 

Lemma 3.3. (Lower bound on the distance between particles, [2D]) 

Let (xi(t))i=i AT he the solution of Ii3.20\) on [0,T] and let 

d{t) := mm{\xi{t) -Xj{t)\,i ^ j} 

be the minimal distance between particles. Then under assumptions Ai) and Aii), we have 
for all t G [0, T] 

d{t) > d{0)e-^' with C = -fK 

where K is the space Lipschitz constant of a. 

This lemma prevents the crossing of particles in finite time. As a consequence, we easily 
get the following long time existence result: 

Corollary 3.4. (Long time existence of a solution to the ODE system) 

Under assumption (A), there exists a unique solution {xi{t))i=i^,,,^N of Ii3.20\) on [0, +00). 
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4 The convergence result: proof of Theorem 11.11 

We want to build a supersolution with a parameter 5 > to fix later. We define {xi{t))i=i^,,,^N 
as the solution of 



-^ = j i-S -a{x^,t) + -y^ —], on (0, +00), 

at \ 71 ^ Xi — Xj I 



x,(0) = x0-5, 

with 7 given in (11. 5p . 
We also define 



for i = l,...,N 



{t,x) = eait,x) + Y.{<f^{^^^)-'~^^(t)^{^^^ 



i=l 

where 

Ci{t) =Xi{t) 

a = ^-^ with a = W"{0) > 
a 

and (j) is given in Theorem 13.11 and tp in Theorem 13.21 

Then we have 

Proposition 4.1. (Supersolution at the initial time) 

Under assumption (A), there exists 60 > such that for all < d < 60, we have 

(4.21) v%0,x)>v'q{x) for all x eR 
for e > small enough. 

Proof of Proposition 14.11 

First we choose 6 > small enough such that 

a;° < Xi+i (0) < x°^i for z = 1, ..., A^ - 1. 
Let A > he large enough such that 

(4.22) -> (X^Iq(O)I) sup IV'I. 

" \1^ J m-A,A] 

Case 1: \x — ^^(O)! > eA for each i = 1, ..., A^ 

/ ^ — ^. ( ) \ f X — X- 
Using the fact that </> is increasing, we deduce that (p I ''- > I 1 and then 



using (14.22p . we get that 

v^{0,x) > Vq{x). 

Case 2: \x — Xi^^{0)\ < eA for some index io E {1,...,A^} 

Then 

X -^io(O) 



e 



> (j)i-A) > 
10 



while 






X -Xio{0) -6 



<(j)iA--\=0(^) as e^O. 



Therefore for e small enough, we have 



> 



X — XI. \ -^r-^ ,_. , / X — Xi{0) 



\ N 



Using again the nionotonicity of 0, we conclude in case 2 that 

v%0,x) >v^q{x). 



Finally case 1 and 2 prove that f l4.2ip holds for any a; G M. This ends the proof of the 
Proposition. D 



Let us define for i = 1 A^ 



^i := 4) 



X - Xi{t) 



X Xi[Z) \ TT I "^ •^iy^J 

- £1 



where H is the Heaviside function. 

Lemma 4.2. (Computation using the ansatz) 

Let 

r ■= evl + - {W\v') - eLv' - ea}. 



Then for any iq G {1, •••, N}, we have 



(4.23) 



r = el + {a~a-a) + 0{ 



"iO 






where t] is defined in Ii3.19\) and the error ef^ is given by 

4 = o{e) + J2 om + E o{4>^) + E ^ 

if^io i^io ij^k) 

Proof of Lemma 14.21 
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We have, using the equation for the corrector, that 



N 



J2-Xi^ 



— ill '^ \ ) 



j=i 



N 



i=l 



^E (^^)'^' 



fix Xiylj 



ecii^i 



e^ot 



e^La — ea 

N 

+W' {ea + J2\<Pi- ^^^'^Pi 

i=l 



N 



i=l 



J2}w'{4>d-ec. 



W"{4>r)i^i + 0' 



fix Xi\Zj 



r^iW"{^;)-W"{Q) 



We first remark that all the terms in (p' vanish (because Cj = ij). Using assumption (A), 
we can bound at, La and q. Moreover using Theorem 13.21 we can also bound ip'. Finally, 
collecting all the terms of order e together, we get simply 



r = 0(e) + 



-ea 

N 

+W' {ea + J2\4>i- eci^i 

i=l 



N 



■EW" 



j=l 



- eci 



W'\^,)xlj, + r^(w'\^i)-W"{S}) 



Now let us choose any index Iq G {1, ..., A^} and let us make a Taylor expansion of W'{ 
when the argument of this function is close to 0jg. We get 



-a + Oie) 



Waj + W"UJ I ea 









-{W'{<j)i^) -eci^ 



W"{4>^o)^,, + v{w"i4>^o)-w"{o) 






w"{^;)iJi + 0(0, 



We remark that the terms in H^'((/)jg) vanish and the terms in ■?/'jo vanish. Now taking into 
account the definition of ef^^, we see that we get 



r = - 



a + e. 



«0 



+ (W^"(0i, 
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and then 



r = el + (aa - a) + 0( 



^^oJ 



a 






VCio 



This ends the proof of the lemma. 

We are now ready to prove the following result: 



D 



Proposition 4.3. (Supersolution for positive time) 

Under assumption (A), there exists 60 > such that for all < d < 60 and given T > 0, we 
have 



(4.24) 



vl>- ( Lr - -W'{v') + a{t, x) ] on (0, T) x 



for e > small enough. 



Proof of Proposition 14.31 

Case 1: \x — Xig{t)\ < e^ for some index zq G {I, •••5^} 

Hence for e small enough, we also have 

\x-Xiit)\>e'^ for te{l,...,N}\{to}. 

Recall also that by ( I3.17p . we have 



4>{y) - H{y) 



We deduce that 



E^ 



aTTy 
1 



< 



C 



for all \y\ > 1. 



i¥=io 



e aiTix — Xi 



< (iV-l)C£3 =0(53). 



Therefore using fl4.23p 



I' > el + 0{e^) + aa - a 



+0( 



I -{6 + a{t, x)) - V — —— 
la ^ an{xi^- 



+ 



--10 



Xi) 07 

> el + 0(£t) + 5 + O(0,J {a(t, x) - a(t, XiJ} 

> s + o{e-^) + el 

> 6/2 for e small enough 

where in the first inequality we have used the fact that a^rj = 1. To get the last line, we 
have used the fact that ef — )■ as e — )■ 0. 
Case 2: \x — Xi{t)\ > e^ for each i = I, ..., A^ 
Using in particular that 

4>i = 0{el) 
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we conclude again from (I4.23P that 

P > 6/2 for e small enough. 
Finally case 1 and 2 prove that for e small enough we have 

r > 6/2 for every x G M 
which ends the proof of the Proposition. 

Proof of Theorem 11.11 

Step 1: existence and uniqueness of the solution 

First remark that for e small enough the initial condition satisfies 

-1<v'q<N+1 
and the functions 

u(t, x) = -1- KH and u{t, x)= N + 1 + KH 
with 

are respectively sub and supersolutions of (ll.3p on [0, +oo) x M. From the Perron's method, 

it follows the existence of a solution v"^ of ( II. 3p on [0, +oo) x M which is moreover unique 

thanks to the comparison principle (Proposition 12. 4p . 

Step 2: sub and supersolutions as e — )■ 

Now, given any fixed T > 0, and thanks to Propositions 14. 1 1 and 14. 3[ we can find some 6 > 

small enough such that W is a supersolution of ( II. 3p on [0, T) x M. 

Similarly we can build a subsolution vf (defined simply as W but with 6 < 0). At the initial 

time we have 

v'{0,x) >v'{0,x) >vf{0,x) for all x eR 

and then from the comparison principle, we get that 

v^t, x) > v%t, x) > vf{t, x) for all t G [0, T), x eR. 

Now for fixed T > 0, we get the result (namely (11.70 and (II. 8p ) on [0, T) x M, passing to the 
limit as e — 7- 0. Finally, because T > was arbitrary, we recover the result for all time on 
[0, +oo) X M. This ends the proof of the Theorem. D 



5 Known and further results on the half-Laplacian 

In this section, we recall several results that will be used in Section 6 for the proof of The- 
orems 13.11 and 13.21 We first recall four definitions (that are equivalent for smooth enough 
functions) of the half-Laplacian operator (Fourier, Levy-Khintchine, classical harmonic ex- 
tension, notion of weak solutions). 

We will essentially use the notion of weak solution (via the harmonic extension) for which 
we will show that it coincides with the notion of viscosity solutions (via the Levy-Khintchine 
formula), when the function is in the space L°°(M) fl Ci^^{M.) D H'H 
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Then we recall results from [8] on the layer solutions and the comparison principle (for 
the harmonic extension). We also recall regularity results for harmonic extensions and give 
new useful results for the L°° regularity of solutions to half-Laplacian type equations. 

There are many points of view to handle half-Laplacian operators. Here we have chosen to 
consider it through a harmonic extension problem, although another more direct approach 
could have been taken. However, our approach has two advantages: first, many of the 
results reduce to its analogous for harmonic functions, and second, we use the setting already 
established by Cabre and Sola Morales in [8]. Another advantage is that this is precisely the 
setting of Cabre-Sire [7] for generalization to any power a G (0, 1). 

5.1 The Direct approaches 

5.1.1 The Fourier approach 

We recall our definition of the Fourier transform for w G S(M.) (where S(M.) is the Schwartz 
space): 

w(^) = j w{x) e-'^'-'dx for all ^ G M 
Jr 
and recall the definition of the following Hilbert space 

H^m) = Le L\R), ^^e (1 + m^m' < +oo| 

with the scalar product 

where v* denotes here the complex conjugate of v. We also set the following Hilbert space 

i7-^(M) = LG5'(M), «;GLL(M), Jj^{1 + \^\)-'\w{0\'<+oo 
with the scalar product 

It is known that iy~2(]R) is the dual of if 2(M). Recall that we set L = — (— A)2 to be the 
half-Laplacian operator defined in Fourier space by 

(5.25) LwiO ■■= -\^\w{0 for ^ eR. 

In particular, L is a continuous linear mapping from if 2(]R) into ii~2(]R) satisfying 



Moreover L is self-adjoint as follows: for any v,w G if2(M), we have 
< Lv,w >^-^(j,),^^(j,)=< Lw,v >^-^(^),^i(^) 
Finally let us remark that for any w G ii2 (M) we have 



(5.26) < —Lw,w > 1 , 1 , , + w"^ = — (w,w) \, ,. 

^ ' ' _ff"^(R)xH2(R) / 27r ^H2(R) 
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5.1.2 The approach in real space 

We have the foUowing resuh (see [28], or the more recent paper |T3], Theorem 1): 

Theorem 5.1. (The Levy-Khintchine formula) 

For anyO < [3 <1, if w e H^R) n C/^f (M), then we have 

Lw{x) = - I —^ iw{x + z) - w{x) - W(x)l|b|<i|) . 

Moreover this formula allows also to define Lw for instance for w E L°°(R) fl C^J^ 

Remark 5.2. In the previous formula the characteristic function l{|2|<i} can be replaced by 
l{\z\<r} for any r > 0, which does not change the value of the integral. 

The operator L is called a Levy operator associated to the Levy measure dz/z"^. In 
particular, the Levy-Khintchine formula allows to see that for any w G H^lM.)^ we have 



(5.27) Hd^ |e||6(0r -< -Lw,w >,.,^,^^,,^,^^llj.dx' 



\x' — xp 



This last formula is a simple exercise and can be found in the book by Lieb-Loss |29| . theorem 
7.12, part ii. 

5.2 The approaches by harmonic extension on M^ and notion of 
weak solutions 

In the following, we give a characterization of the half-Laplacian operator in R through 
a harmonic extension to the upper half-plane as a Dirichlet-to-Neumann operator. This 
relation was well known but it has been recently rediscovered by Caffarelli-Silvestre in [9], 
where they also give the generalization to any fractional power of the Laplacian between 
zero and one. 

The main advantage of this interpretation is that the regularity results and maximum 
principles for the half-Laplacian can be shown through the analogous results for harmonic 
functions in the extension, as it has been considered by Cabre and Sola Morales in J8]. 

5.2.1 Classical harmonic extension for smooth functions with compact support 

Assume that w G C^(M) (the space of smooth functions with compact support) and consider 
its harmonic extension u{x, y) on 

Q:={{x,y): x G M, y > 0} = M^ 

defined as the solution of 

{Au = on Q, 

u{x, 0) = w{x) for all x G dVl. 

Then it is known that 

du 
(5.29) Lm;(x) = - — (x,0) for ah x e dVt. 
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Here, and for the rest of the paper, v denotes the exterior normal along 90. 

The half-Laplacian appears to be a Dirichlet-to-Neumann operator. This is also called 
a Steklov-Poincare operator. Moreover, we can use the Poisson kernel P for the upper 
half-plane to write u = P *xW, i.e., 

(5.30) u{x,y) = C2 -2 -w{Od^, 

Jk \x - Q + y' 

where C2 is a dimensional constant. 



5.2.2 Harmonic extension of functions in H2 

Now let us define the scalar product 



{u,v)h = / Vu ■ Vw + / uv 
Jn Jan 

and consider the Hilbert space 

Then we have the following result: 

Theorem 5.3. (The trace mapping) 

The trace mapping 

T: H — > H^R) 

u I — )• Tu := u{-, 0) 

is linear continuous and satisfies 



(5.31) \\Tuf 1 <27r||n 



|2 

\h ■ 



Moreover we have \\Tu\\ 1 = 27r ||n|| r. zf and only if u is harmonic on Q. 
II iih3(r) II ii-H ■' *" ■' 

Proof. For u smooth and with compact support, the classical embedding trace inequality 
reads 

f , ,^, iTuj x') - Tu ix)\'^ 
\x' — x\ 



(5.32) ||Tm|| . 1 := / dxdx' — < 27r / IVmI dxdy. 

ri^(M) /.^ „ \nr*' — nrl'^ 



A proof of the previous inequality can be found in corollary 6.3 in j2]. Then (I5.3ip follows by 
putting together the previous line with ( 15.26P and ( 15.27p . For a general u the result follows 
by density, while the equality case is dealt with in Theorem 15.41 below. 

Now, we can reformulate problem (I5.28p -( l5.29p in a more general setting: 

Theorem 5.4. (Half-Laplacian defined by the harmonic extension) 

For any w G H2(M.), there exists a unique u & H which is harmonic on Vt, such that Tu = w, 
and it is written as u = P *xW. Moreover for any v G H , we have 
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In particular, 

and 

(5.34) I Nuf + w^ = — \\wfi^. 
^ ' Jn ^ Jan Stt " "^Mm) 

The previous theorem allows to define weak solutions (in the sense of distributions) for 
half-Laplacian equations, that may not be necessarily equivalent to the notion of viscosity 
solutions introduced in Definition 12. 1[ Although the introduction of a new concept may seem 
confusing at first, it will allow to use the available functional analysis tools. Moreover, this 
definition is a particular case of the one introduced in [8]. Thus we set: 

Definition 5.5. (Weak sub/super/solution) 

Let d be a measurable bounded function on M, and let g G L^(M). We say that w G H^^R) 
is a weak supersolution for the equation 

(5.35) Lw — d{x)w = g on dfl 
if, for u := P *x w and 



B{u,v):= / Vu-Vv+ / d{x){Tu){Tv), 
iQ. Jan. 



we have that 



(5.36) B{u,v)> gTv forallvEH. 

Jan 

Respectively, w is a weak subsolution if 

(5.37) B{u,v) < I gTv forallveH. 

Jan 

Accordingly, we say that w is a weak solution if it is both a weak sub and supersolution. 

We summarize the different concepts in the following lemma: 

Lemma 5.6. (Equivalence of solutions) 

Given w G L°°(]R) fl Cj (M), the following two ways of compuing Lw, i.e, 

1. as a Dirichlet to Neumann operator (see section 5.3.1., in particular (I5.28p - (l5.29p ).• 

2. using the Levy- Khintchine formula (see theorem \5 . 1\) . 

are equivalent. 

If, in addition, w G H^^R), then the half-Laplacian definition in term of distributions as 
given in theorem \5.4\ is also equivalent. As a consequence, the notions of viscosity and weak 
solutions from Definitions lKj\ and \5.5[ respectively, coincide in the smooth case. 

Proof. We refer the reader to the paper |13] . 



5.3 Miscellanea 

Here we remind the reader some basic facts. 

For harmonic extensions, we have (cf. |35], chapter II. 2): 

Theorem 5.7 (L^ bound). If w E ^^(IR), 1 < p < C)0, and u = P *x w is the Poisson 
integral of w, then u is harmonic in M^, Ymiy^ou^x^y) = w{x) for almost every x eM and 

M-^y)\\LP{R) = I Hx^y)]" dxj < ||w||ip(K) 
for all y > 0. If w is bounded on M, then also u is uniformly hounded in M^. 

As a consequence, given w G -ffz (M), then for u = P*xW E H we have that u G W^''^{B^) 
for every half ball B'^ , since 

u^ dxdy < / u^ dxdy < R w"^. 

B+ Jye{0,R) JxeR JdQ 

Moreover (cf. [T7], Chapter 2, Section 2, Theorem 7): 

Theorem 5.8 (Local gradient estimates). Assume that u is harmonic in a domain U. 

Then 

C 

|/^"M(xo)|<^||n||^^(5^(^^)), 

for every ball B^.{xq) C U , and multiindex \a\ = k. 

Next, we give some classical facts on Sobolev spaces, from [30] (Section 7.1.3, Theorem 
4, part (6) and Section 8.8.1, Theorem 4): 

Theorem 5.9. (Sobolev embedding) The injection 

is continuous, for any p > 1. Moreover, if K G R is a compact subset, then the restriction 
map 

H^(R) ^ LP{K) 

is a compact embedding for any p > 2. 

5.4 Layer solutions for half-Laplacian problems 

The paper [8] is concerned with the boundary reaction problem 

{Au = in Q, 

— = f{u) on dn, 

for some non-linearity /. We say that m is a layer solution of (I5.38P if it satisfies (15. 38 p . 
Ux > on dVt and 

lim u(x,0) = 0, lim m(x,0) = 1. 

a;— >■— CO a:— >+oo 
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Since we are only concerned with smooth solutions of ( I5.38p . then we do not need to 
worry here about the different concepts of weak or viscosity solutions for half-Laplacian 
problems previously introduced. 

In the next theorem we summarize the basic facts about existence and behavior of layer 
solutions for f l5.38p (see theorems 1.2, 1.4, 1.6, and equation (2.27) from j^). 

Theorem 5.10. (Known existence and properties of the layer solution) 

Let f be any C^'^ function with < (3 < 1 and such that W = — /. Then there exists a layer 
solution uw G C, (fi) of (I5.38P if and only if 

W'{Q) = W\l) = Q and W >W{'d) = W{l) zn(0,l). 
Moreover, 

1. // Vr"(0), Vr"(l) > 0, then a layer solution of (I5.38P is unique up to translations in 
the X variable, and it satisfies 

(5.39) {uw)x > inTl. 

2. We have the estimate 

(5.40) \Vuw{x,y)\<-—-. -, forallxeR,y>0. 

1 + \{x,y)\ 

3. Set (j){x) := uw{,x,Q). We have the bounds 

C 

(5.41) \(j)(x)-H(x)\< — forallxeR, 

1 + |x| 

C C 

(5.42) ^ < 0'(x) < ^, for all x G R, 

1 + |x| 1 + |a;| 

for some positive constants C, Ci, C2. 
4- A layer solution is always a stable solution. This means that 



(5.43) f \Vvf + f W'icpy > 

Jn Jan 

for every function v G C^{Q) with compact support in Vt. 



One of the crucial ingredients in the present work is the precise knowledge of the asymp- 
totics for as given in (I3.17p . In order to improve from the known ( 15.4ip . we will use a 
comparison argument with the explicit layer solution of the Peierls-Nabarro problem (see j8]. 
Lemma 2.1, for the proofs, or the original paper [36]): for each a > 0, consider the equation 

An" =0 in (], 

-Q^ = -K(w ) on on, 
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with potential given by 

K(t) = -^sin[2vr(t-l)]. 

It is well known that the function 

(5.45) m"(x, y) = — arctan h - 

71 y + a 2 

is a layer solution for the problem fl5.44p and, in fact, it is the only one. 
A simple computation gives that 

(5.46) < ^ «^ + 2< 



()"(x) - H(x) + — 

TXX 



< 7T for all X G 

1 I I I ^ 

1 + \x\ 



along dQ. Set (j)°'(x) := u^{x,0). We know explicitly its asymptotic behavior 

(5.47) 

an improvement from (I5.4ip . 

Finally, fix u := {u\y)x E C, (fi), and note that u = P *x <P' ■ Then: 
Lemma 5.11. (Properties of 0'). For u defined as u = P *^. cf)' , we have that 

1. ueH. 

2. (f)' is a weak solution of Lw — W"{(f))w = in the sense of Definition \5. 5i i.e., 
(5.48) B{u, v) = for all v G H. 

Proof. We recall estimate (I3.16p . It shows that (p' G L^(]R). To complete the proof of 
the first claim we need to check that J^ |Vn| < oo. 
Step 1: I{o<y<R} IV^P < +00. 
We proceed as in the proof of Lemma 2.3 in [8]. We first set 













v{x,y) = / m(x, 
Jo 


t)dt. 




We 


can 


see 


that V 


is 


a (smooth) solution of 

r Av = W"{4>{x))(P'{x) 
1 v{x,0) = 


in 
on 


n, 
on. 



We know (I5.40p which implies that v E L'^ {M. x [0, R]). Because the right hand side of the 
equation is bounded in L^(]R x [0, R]), and here we are using again the precise decay of </)', 
then we get also that v G VF^J^ using the boundary gradient estimates for harmonic functions 
in [23], Lemma 9.12. As a consequence, u G iy^'^(M x [0,-R]). 
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step 2: j^y^j^^ |Vm|2 < +00. 

Far away from the boundary of f], we can use standard interior gradient estimates (Lemma 

I5.8J) for the harmonic function -u, say, at a ball Bi{x,y) G Q so that 

where we have used (I5.40p in the last inequality. 
Step 3: Conclusion. 

From Steps 1 and 2, we get that u E H . 
Step 4: Proof of (54^. 
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From fl5.38p with / = -ly', we get for any w G C^{Q): 

= I Vu-Vw+ I wW'{u) 
Jn Jan 

Setting w = Vx with v G C^{Q), we get by integration by parts 

B{u,v) = 
Finally, we conclude by density of C^i^Q) in H. 

5.5 Comparison principles for smooth enough solutions 

Next, we consider some comparison results for half-Laplacian equations: 

Lemma 5.12 (Strong comparison principle. Lemma 2.8 in 0). Let v be a bounded 
function in fl, v E C^(r2) and C^ up to the boundary dQ, satisfying 

{- At; > inQ, 

dv 
— + d(x)v > on on, 
ou 

for some bounded function d satisfying d{x) > r for all x gM. and some r > 0. Then f > 
in Vt unless f = 0. 

We will also need the following variation: 

Corollary 5.13. Let v be a bounded function in Q, v E C'^{Q) and C^ up to the boundary 
dQ, satisfying 

Av = in Q, 

\ dv 

- — h d{x)v = M{x) on dQ, 

where d is a bounded function satisfying d{x) > r for a// x G M and some r > 0. Assume, 
in addition, that the right hand side M can be estimated by 

C 

|M(a;)| < 2' ^^^ 

1 + \x\ 



22 



for some constant C > 0. Then there exists a non-negative constant C depending on C,t 
such that 

c 



f X, 



< 



\x\ 



for all X G 



Proof. Fix a := 2/r, and consider the model layer solution m" given in f l5.45p . with trace 
(j)"-. Because of hypothesis on M and the decay of {4>°')x from (I5.42p . we can find a positive 
constant K such that |M(x)| < K(j)'^{x) for all x G M. Next, compute 



dv 
du 



+ d{x)v = M(x) < K<. < K 



-— ^ + 2< 



< 



2K 



T 



du 



d{x)ul 



along 9n, 



where we have used ( I5.46P in the second inequality. Now we apply Lemma 15.121 to the 
difference [^m^ — v\ to get v < ^m^ in Vt. The same argument, applied to —v gives 
that \v\ < — m". Taking into account the estimate (I5.42p . we achieve the conclusion of the 
corollary. 

D 



5.6 Further properties of weak solutions 

The aim of this subsection is to extend the classical results for harmonic functions to the 
boundary Sfi, in order to obtain regularity and maximum principles for weak solutions (in 
the sense of definition 15. Sp . We set 

ti~ := — infJM, 0}, m"^ := supJM, 0}. 

Our first result is a weak maximum principle for non-smooth solutions: 

Proposition 5.14. (Weak maximum principle). Let u G H be a weak supersolution of 
fl5.35p with g = 0, and let d{x) be a measurable bounded function. Set u = P *xW. Assume 
that u > in B^, and that d > on M\[—R, R]. Then 



u > in 



^+' 



and also, Tu > on 



Proof. This is a small variation of the classical proof for the Laplacian, that can be 
found in Gilbarg-Trudinger [23], Theorem 8.1. Note that {Tu)~ = T{u~) a.e. in M. Set 

VL~ := supp{u~) and 6 := supp{{Tu)~) C M. 

We know that Q~ is contained in M^\i?^ and O C M.\{—R, R). Use u~ E H as an admissible 
test function in f l5.36p : then 



< - 



|Vnr- / d{x)u' 
'e 



The proof is easily completed. 
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D 
Note that a function u & H must have trace Tu belonging to Lp(M) for all p > 1. 
However, it may not belong to L°°(]R). In the next result, we use a Moser iteration scheme 
to show that weak solutions are actually bounded. 

Proposition 5.15. {L°° bound for g = 0). Let w E H^iW) be a weak suhsolution of fl5.35p 
with g = for some uniformly bounded function d defined on M. Then 

If w is a supers olution, the conclusion reads 

sup{-u;} <C\\w^\\^^^^y 

The constant does not depend on w. As a corollary, u := P *x w is also uniformly hounded 
on VL. 

Proof. This is the classical proof by Moser for the Laplacian case. We follow closely the 
arguments from Theorems 8.15 and 8.18 in Gilbarg-Trudinger [23], so we will not give every 
single detail here. Assume that w is a subsolution, and set u := P *x w. Given cr > 1, we 
insert into (I5.37p . the test function 

V = {u+y. 

Note that Tv = {Tu)" a.e.. Here we point out that this v might not belong to H\, in order 
to make the argument rigorous, a truncated version must be used instead in the definition 
of w, say replacing u^ by 

(5.50) minju'*', ra}, for n G N. 

We obtain 

a I u"-^Vu-Vu^ < I -d{x){u^Yu. 
Jn Jdn 

Next, we set Ui := {u~^y and Wi := Tui for 7 = ^y^, and use that \d{x)\ < Cq. Thus we get 



/ iVuil^ <Ci(7) f wl with Ci 
Jn Jdn 



(5.51) / iVuil^ <Ci(7) / wl with Ci(7) 



27-1 



For the left hand side, we proceed as follows. We use the trace inequality (15.321) to 
estimate 

■'■ II l|2 ^ /" IV7 |2 

\\Wi\\ .1 < / VMi . 



But, on the other hand, by Theorem 15. 9[ for every p > 1, there exists a constant Cp > 
such that 

<C2(7)lki|li2 with C2(7) = Cp(l + 27rCi(7)), 

where we have used (I5.5ip in the last line. 

If we denote 

1 

^(g):= I / \w+\'^' 
\Jdn 
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then we have shown that 

^(2p7)<[C2(7)]^*(27)- 
Fixed p > 1 and choose 7^ = p™ for ?Ti G N; by iteration we know 



m—l 

(5.52) ^(27^) < C„*(2), VmeN\{0}, C„ := J] [^'2(7^)]'^ • 

j=0 

It is then easy to check that Cm < Coo uniformly bounded in m G N (which is also shown, 
for instance, in page 190 of j23])- Remark that 

\ +\ r f ^ f I +1"^' 
sup \w \ = hm I — / ky^ 

{-R,R) 1^°°\'2R J{-R,R) J 

Letting m to infinity in (I5.52p . we have 

sup \w^\ < Coo ||^i'+|L™^> 

for every R> 0. We get the first part of the conclusion, taking R — > +00 and then 7i — )■ +00 
in (I5.50p . The second part is analogous, using u~ instead. The final assertion is an immediate 
consequence of Theorem 15.71 

D 
We will use also the following particular case (that is not the optimal one can prove, but 
it suffices for our purposes): 

Corollary 5.16. {L°° bound when g is bounded). Let w G H2(R) be a weak solution 
of (I5.35P with non-vanishing right hand side g. If there exists po > 1 such that 

(5.53) ||5'IIlp{r) < c for all p > po 
(with c independent on p), then w is bounded. 

Proof. We use Proposition 15.151 with some modifications in order to account for the 
extra right hand side term involving g. By hypothesis, u; is a weak solution satisfying for all 
veH, 

B{u, v) = — gv, ioi u = P *x w. 
Jan 

Similarly as before (using if necessary a truncated version), we use the test function v = 
\u\'^~ u, Ml = |m|'^~^m and Wi = Tui. We use Holder inequality with exponents 27, y^^ 

which gives 

1 ^ ^ 27-1 



gv 
an 



we get 



an J \Jdn 



2^"" .... r^UA-^!^.f f ini^^^'-' 



|V«ir<Ci(7)/ wi + C,{^)[ wt] with Ci(7) = :-^- / \9 

Jan \Jan J 27 — 1 \Jq^ 
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Assuming (I5.53P with 2j > pq, we get 



iVnil" < Ci(7) max <{ 1, / w'f)- with ^1(7) = ^1(7) + ^^ 



'n I Jan J 27 - 1 

The rest of the argument follows similarly, with 6*2(7) = max |l, Cp(l + 27rC'i(7))}, 

^{q) = maxjl, (/g^ |u;|'')«|, and ^(2^7) < {C2{'y))^^{2'y), starting the iteration at 

27 = Po- □ 

We now quote the regularity theory for bounded weak solutions of [8]. 

Proposition 5.17. (Regularity of weak solutions). Let R > 0, (3 E (0,1), and let 
w e H^R) be a weak solution of IK35\f with d G C^([-4i?,4i?]), g G C^([-4i?,4i?]). // 
u = P*^w E L^iB+j^) n iyi'2(5+j), then u e C^'^{B%) and 

for some constant Cr that depends only on P, R and on the upper bounds for \\u\\^oo(q+ ), 

IMIIc/3([-4R,4R]) ^'^^ llfi'llc''([-4i?,4i?J)- 

Proof. This is a consequence of Lemma 2.3 from j8], using definition 15.51 for weak 
solutions. 

D 
In particular, we obtain C, regularity up to the boundary dQ: 

Corollary 5.18. (Uniform regularity of weak solutions). Given R, (5, w, u, d, g as 

above, assume, in addition, that \\d\\^fi,^^ < C, \\g\\ci3(M.) ^ C. Then u G C^'^iM. x [0,i?]) and 

\\^'^\\l°°{Rx[o,r]) — ^Ri 
for some constant Cr that depends only on 13, R and on the upper bounds for \\u\\jj, ||rf||c/3dR) 

(^nd ll^llc^(K). 

Proof. We apply the previous proposition to every ball Bji{xQ, 0), Xq G M, as in Lemma 
2.3, part (b) of [8]. It uses the condition that u is uniformly bounded on ^2, which in our 
case follows from Corollary 15.161 

D 

6 Results on the layer solution and the corrector 

In this last section, we give the proofs of Theorems 13.11 and 13.21 

6.1 Proof of Theorem 13JJ 

Let W he a potential satisfying hypothesis Ai). Theorem 15.101 shows existence, regularity, 
uniqueness and asymptotic behavior of a layer solution < Uw < 1 satisfying Mvy(0, 0) = 1/2, 
u E C-,'^ (fi), for the problem 

An = in r2, 

7— = -W (u) on dQ. 
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We denote its trace value by := Uw{-, 0). The aim of the next paragraphs is to show the 
more precise estimate (I3.17P for the asymptotic behavior near infinity. 

Set a = l/a, where a = W"{0) as defined in fll.6p . and choose the model layer solution 
u"' as given in fl5.45p . We know explicitly because of f l5.47p that 0" := u"(-, 0) has the right 
asymptotic behavior (I3.17p . We would like to compare u^ to u"", in order to obtain the 
analogous estimate for uyy. Therefore, we set v := Uw — u"-. Then v satisfies 

Ay = in fi, 

^ = -W'{<P) + V;{r) on dQ. 
Now we estimate 

-W'{<P) + W\r) = -W"iO)[<p -01+0 (-^) , 

VI + \x\ J 

thanks to the growth estimate (15.4ip . And also, 

-w'ir) + K'(r) = [-w"{o) + v^'mr + o (—^) . 

\l + |x| / 

Because of the choice of a, we have that iy"(0) = V"^"(0). Thus adding the previous two lines 
we conclude that f is a solution of 

Ai; = in fi, 

dv 

— + W"{<d)v = M{x) on dn, 

for some function M satisfying 

VI + \x\ 

Therefore, fl3.17p follows from Corollary 15.131 (comparison principle for smooth solutions) 
and the fact that W"{Q) > 0. 

The proof of Theorem 13.11 is completed due to the fact that a smooth, bounded solution 
u of f l6.54p . with trace 0, is a viscosity solution in the sense of Definition 12. II of 

L0 - ly' (0) =0 in M, 

with L defined through the Levy-Khintchine formula, as it is noted in Lemma 15.61 D 

6.2 Proof of Theorem 13.21 - coercivity 

Here we remind the reader of the notation introduced in Section |5l let be the layer solution 
constructed Theorem 13. H with harmonic extension Uw (some properties are given also in 
Theorem 15. lOp . Its asymptotic behavior is given in Theorem 13.11 



(6.55) 




x) — H(x) 



C 
< 



anx 
|0'(x)| < 
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1 + x^' 
C 
l + a:2' 



for all X G M. We also fix u := {uw)x ^ C, (fi), and note that u = P *x <P' ^ H . For some 
of its properties, see also Lemma 15.111 



In the following, we set up the necessary functional analysis ingredients for Theorem I3.2[ 
while the proof will be completed in the last subsection. 

We set Hq to be the subspace in H consisting of functions whose trace is orthogonal to 
(j)' in L'^idn), i.e. 



(6.56) Ho:= <ueH : u(j)' = 0}. 

I Jdn 

Lemma 6.1. Hq is a closed subspace of H . 



Proof. As we have noted in Theorem 15.31 the trace of functions in iiT is a well defined 
function in H2(M.). In particular, the trace mapping T : H ^ L^(]R) is continuous by 
Sobolev embedding (Theorem 15. 9p . Then, if we have a sequence {uk}'^i -^ Uq in H with 
Uk & Hq, then 



an 



1/2 / /• \ 1/2 

'an J \Jan 



< ( / 0'^ I I / \uk — Uo\ ] — 7- when /c — )■ oo. 



using the bounds for 0' given in (I6.55p . We conclude that Uq G Hq, and the lemma is shown. 

D 
Now, given u,v & H, we define the bilinear functional 

(6.57) B{u,v):= Vu-Vv+ W"{(j))uv, 

Jn Jan 

whose associated quadratic form is precisely 

Q(u) = f |Vn|2 + f W"{(j))u\ ueH. 
Jn Jan 

It is easy to check that 5 is a continuous and symmetric bilinear functional in H. 

The key lemma in this section is the following: 

Lemma 6.2. (Lower bound for Q). There is some constant Co > such that for all 
u G Hq, where Hq is defined in (I6.56p . we have that 

Q{u) >Cq [ u\ 
Jan 

So we immediately have that 
Corollary 6.3. (Coercivity). There exists some C > such that for all u G Hq, 

B{u, u) > C \\u\\^ . 
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Proof. It easily follows from Lemma [6.21 Indeed, 

B{u,u)= [\Vuf+ f W"{(l))u^>{l-^i)Co I u^ + f^Q{ 
Jq Jdn Jan 



> /i / I Vnl' + / u' [(1 - i2)Co + /iW^"(0)] . 
Jn Jdn 

Because W" is bounded, there exists // > small enough such that (1 — fi)Co + iJ.W"{(f)) > ^. 
Therefore, 

(6.58) B{u,u)>mm{fi,^}\\u\\%, 

as claimed. 

D 

Proof of Lemma 16. 2t It is divided into several steps: 

Step 1: Stability. Because of the construction of uy/ and </>, we know from Theorem 15.101 
that Q{u) > for all u G C^{Q) (i.e. layer solutions are stable). Thus the functional Q is 
bounded from below also in H. 

Step 2: Existence. 
Step 2.1: Preliminaries. 
We seek a minimizer for the functional 

Q(u) = [ \Vuf + / W"i(P)u'' 
Jn Jdn 

subject to the constraints 

u E Hq and / n^ = 1. 
Jdn 

Note that W"{(j)) is not bounded from below by any positive constant. However, due to our 
assumptions Ai) on W and the growth condition (I6.55P for </>, we still can find some r > 
and h a non-negative smooth function supported on a compact set K := [—R, i?] C M that 
makes W"((j)) + /i > r > for all a; G M. We write 



Qiu) 



[iVuf+f iW"i(l)) + h)u^ - f hu" =■ Qi{u) + Q2{u 
Jn Jan J Jk 



Let {tifcjfcii C -f^o be a minimizing sequence for Q with Jqqu\ = 1. This L^ bound for u^ 
tells us that, in particular, \Q2{uk)\ < sup /i < C for all k, for some constant C > 0. As 
a consequence, we also know that there exists a constant C such that Qi{uk) < C . From 
here we get a uniform bound for the if-norm of the sequence {ttfc}, depending on r and C . 

Next, we concentrate on the space ilf2(M). We define an equivalent Hilbert space norm 
||-||^ on M as follows: for any w G iJ2 (M), we consider its harmonic extension u := P *^. w, 
and define 

\\M\h ■= Qii^)- 
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We have shown that ||Mfc||^ < C for all k. This implies that for Wk '■= Tuk, the norm 
ll'U^fell// < C' for another constant C. Because if is a Hilbert space, then there is a subse- 
quence (still denoted by Wk), such that Wk -^ Wq, weakly in H, for some Wq G H2(M.). Weak 
convergence implies, for instance ({5|, proposition III. 12), that 

\\wo\\f[ < liminf ||u;fc||H- 

fc— >oo 

Set Uk := P *x Wk, that has the same trace as Uk, and uq := P *x Wq. It is clear that 

(6.59) Qi{u,) < Qi{u,) 

and that equality holds when Uk is already harmonic. Then we have shown that 

(6.60) Qi{uo) = \\wo\\jj < liminf ||wfc||^ = Qi{uk) < Qi{uk). 

fc— >oo 

On the other hand, the embedding if 2(]R) --^ L^^K) is compact for any compact subset 
K (Z M. and p > 2, as stated in Theorem 15.91 Thus, we can find some subsequence (still 
denoted by Wk) such that Wk — ?■ Wq strongly in U\K). As a consequence, also Wk — )■ Wq in 
L'^[K). The uniqueness of the limit will give that Wo = wq, at least a.e. in K. In addition, 
since Q2 lives only on the compact subset K G M., we also have that 

(6.61) Q2iuk)-^Q2iuo). 
From (I6.60p and (I6.6ip we deduce that uq G H satisfies 

Q{uo) < \imini Q{uk). 

fc— >oo 

Remark that if lirainik^ooQiuk) > 0, then we have already finished the proof, because 
this implies (I6.58p . Therefore we can assume that 

(6.62) liminfQ(ufc) =0 

Step 2.2: Concluding. 
Case A: Uq = 
Remark first that 



Mfc = r > 0. 
Ian 



Qi{uk) >r 
On the other hand, we have Q2{uk) — )■ Q2{uq) = 0. Therefore we get 

limini Q{uk) = liminf {(5i(nfc) + Q2{uk)} > r > 0, 

which is in contradiction with (16.621) . Therefore only the next case occurs. 

Case B: uq ^ 

Then we have < Qi{uq) = and we will show that Tuj. — )■ Tuq in iJ2(]R). Since {uk} was 

a minimizing sequence, then Q{uk) — )■ = Q{uo) when k — )■ 00. In addition, we have shown 

that Tuk — 7- Tuq strongly in L'^{K), so that Q2{uk) -^ Q2{uo). Because Q{uq) = < Q{uk) < 

Qi^k) — > 0, we deduce that Q{uk) — > Q{uq). Moreover we have Q2{uk) = Q2{uk) -^ ^2(^0) 

and Q = Qi + Q2 which implies that 

(6.63) Qi{uk) -^ Qiiuo). 
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By definition (I6.63P means that 

(6.64) lim \\wk\\H = Ikollj^- 

fc— >oo 

Proposition III. 30 in |5] assures that if we have a weakly convergent sequence Wk —r Wq in 
the space H and such that we have convergence of the norms fl6.64p . then the convergence 
Wk — )■ Wo is strong in H, and then in H^[M.). 

As a consequence, Wk = Tuk — )■ Tuq = wq in L^(]R), which, in particular, implies the 
non-degeneracy of the L'^{dQ) norm of wq: 

wl = lim / wl = l. 

an ''-^°° Jan 

To finish, just remark that strong convergence in L^(M) implies that Jq^UqCJ)' = 0. Thus, 
Mo e Hq. 

Step 3: Euler- Lagrange equation. Once the minimizer mq is found, the theory of Lagrange 
multipliers for minimization problems (see, for instance, [T7], chapter 8.4) gives that Uq & H 
must be a weak solution of 

(6.65) B{uo,v) = A / uqv + /x / (l)'v for every v G H, 

Jan Jan 

where A, /i G M. The value of the multipliers can be found by choosing the right test function 
V. Thus if we test the equality (16.651) with v = P *x 4>', we obtain /i = (see Lemma [5. lip , 
while testing against v = uq gives that 

(6.66) A = Q(no)= [ \Vuof + [ W"{<f))ul (> 0). 

Jn Jan 

Step 4: Regularity. First, Proposition 15. 15l with d{x) = W"{(f){x)) — A implies a uniform 
L°° bound for both wq and uq = P *x Wq- Then, Proposition 15.171 assures that uq E C, (fi), 
with the following boundary estimate 



<c(/3,i?,holLo.(^|^),||W^" 



ll"0|lci.''(B+) - ^ [P^-^^ II^0|Iloo(s+^) , \\VV [(P)\\c0{[-4RAR]) 

Note that mq is as smooth as we want in the interior. 

Step 5: Positivity. In the following we show that if A = 0, then Wo = Tuq must necessarily 
be a multiple of </>', which is a contradiction with the fact that uq G Hq and Jq^ Uq = 1. 

Thus we take A = in (I6.65p . Then we have a solution Uq satisfying B{uo, v) = for all 
V E H . We remind the reader that have defined Uy/ to be the layer solution for our potential 
W , with trace := Uw{-i 0), and u := d^uw-, which is also uniformly bounded and has trace 
'u(-,0) = 0'. Because of the hypothesis on the potential W and the asymptotic behavior of 
the layer solution 0, there exists some r > such that W"{(j)) > r on ]R\[— i?, i?] for some 
i?>0. _ 

We have shown, in particular, that Uq is uniformly bounded in Vt. Since u is positive 
from (I5.39p . then there exists k > such that |uo| < t^u in B^. We set v := nu — Uq, and 
note that v satisfies the conditions of Proposition 15.141 Thus we can conclude that i; > in 
M.\ so that Uq < Ku. An analogous argument gives that —Uq < ku. We have proved that 

31 



|mo| < nu on R^ and, in particular, |w(x)| < k0'(x) for all x G M. 

The next step is to decrease k up to the first point where ku touches |mo|: assume, without 
loss of generality, that for some k* > there exists (xo,yo) ^ ^\ such that 

(6.67) k*u = Uq at (xo,|/o) and k*u>Uq on Vt. 

We set di := max{iy (0), r}, and ^2 := di — W"{(j)). We note that both are non-negative 
functions on M and di{x) > r for all x G M. Let v := k*u — Uq > 0, that is a weak 
supersolution of 

At; = mil, 

dv 

— — h div = d2V > on dQ. 

ou 

We are in the hypothesis of Lemma I5.12[ that can be applied because the functions have 

enough regularity (in particular v G C, (fi)). However, because of (I6.67P we necessarily 

must have t> = 0. We have shown that Uq = n*u, which at the boundary translates as 



w = K*(h'. 



Step 6: Conclusion. The multiplier A in (I6.66P must be non-negative. However, we have 
shown that A 7^ 0, so then 

(6.68) Q{u) > Q{uo) = A =: Co > for all u e Hq with / n^ = L 

Jan 

Finally, substituting 

u 
Ml := 



'Ian-' 
for any u & Hq into (I6.68P finishes the proof of Lemma 16.21 

6.3 Proof of Theorem 13.21 - regularity 

Here we seek a solution for 

(6.69) L-^ - W"{(f))ilj = g 

where 

g:=<P' + ^{W"i<P)-W"m with V = ^^y 

It is easy to check that 

geH^R), [ (70' = 0. 
Jan 

We define also the continuous linear functional on H given by 



D 



Ev := — gTv. 
Jan 

The previous Corollary 16.31 asserts that the bilinear functional B from (I6.57P is coercive on 
the space Hq. Then, Lax-Milgram theorem can be applied to obtain a unique weak solution 
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Uq ^ H such that B{uq,v) = Ev for all v & H. Moreover, the solution depends continuously 
on the initial datum: 

IIMoII rr < C llgll 1 . 



We set ip ■.= Tuq G H^iW). 

The following theorem summarizes its regularity and properties: 

Theorem 6.4. The solution ^ E H^R) of / fOPI) has regularity Cj^'^(M) H L°°(]R), satisfies 
^(±cx)) = and ||V''IIl°°(r) < °^- 

Proof. Note that uq E H so, in particular, ?/; G i/2 (M) (and is also in all L^(M) for p > 1). 

Claim 1: ^ is bounded; as a consequence, also uq. This follows from Corollary 15. 161 after 
estimating the norm of g. We use the asymptotic behavior and boundedness of 0, 0' from 
Theorem I3.H and the Holder regularity of W" . Then 



r< / c^+ / ^<c^ 



2i? + 



(p/3 - 1)RpP-^ 



{-R,R) Js.\i-R,R) ixf'^ 

In particular, for i? = 1, and p > ^ we have, 

II5'IIlp(r) — -^'^• 

Claim 2: C, (M) regularity for t/j follows from Proposition 15.17^ while the uniform 
bound for tp' is a consequence of Corollary 15.181 For the values of ip at infinity, just note 
that V' is a smooth function in L^(M). 

D 

Finally, a weak solution ip e H^R) D C^^^{R) n L~(R) of iK69\i is actually a viscosity 
solution of the equation ( 13.1 8p . thus completing the proof of Theorem 13.21 This is so because 
of Lemma 15.61 
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